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Abstract. A GBDT version of the Backlund-Darboux transformation is constructed for 
■ a non-isospectral canonical system, which plays essential role in the theory of random matrix 

models. The corresponding Riemann-Hilbert problem is treated and some explicit formulas 
are obtained. A related inverse problem is formulated and solved. 
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^ : 1 Introduction 

O . 

t^- , We shall consider a non-isospectral system 

O 

w x (x, z) = i\JH(x)w(x, z), A = (z — x)^ 1 , (1-1) 

^Cj ■ where w x = -jzW, and J and H(x) are m x m matrices: 

03 ■ 

H(x) = H( x y, j = r = j- 1 . 

When the Hamiltonian H > 0, and the spectral parameter A does not depend on x, the system 
^ ■ above is a classical canonical system. A version of the Backlund-Darboux transformation (BDT) 

for the classical canonical system have been constructed in [15]. In our case (jl.lh the spectral 
parameter A = (z — x) -1 depends on x, and here we construct BDT for this case. 

BDT is a fruitful approach to obtain solutions of the linear differential equations and systems. 
It is also widely used to construct explicit solutions of integrable nonlinear systems. For that 
purpose BDT is applied simultaneously to two auxiliary linear systems of the integrable one. 
BDT is closely related to the symmetry properties. Since the original works of Backlund and 
Darboux, a much deeper understanding of this transformation has been achieved and various 
interesting versions of the Backlund-Darboux transformation have been introduced (see, for 
instance, [H [2l [6l [TT1 [T2| [T3l [2T\ I23j). Important works on the Backlund-Darboux transfor- 
mation both in the continuous and discrete cases have been written by V.B. Kusnetzov and his 
coauthors (see [91 [10] and references therein) . 

We apply BDT to construct explicitly new solutions of the Riemann-Hilbert problem on the 
interval [0, I]: 

W + {s) = W^{s)R(s) 2 , 0<s<a, (1.2) 



*This paper is a contribution to the Vadim Kuznetsov Memorial Issue 'Integrable Systems and Related Topics'. 
The full collection is available at |http://www.emis.de/journals/SIGMA/kuznetsov.html| 
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where W(z) is analytic for z ^ [0, a], and W(z) — > I m , when z — > oo, I m is the m x m identity 
matrix. For important classes of R the solution of problem (|1.2|) takes the form 

W(z) = w(l,z), W + (s) = lim w(l,s + irj), W_(s) = lim w(l,s — irj), (1.3) 

77— >+0 r?— >+0 

where the m x m fundamental solution w of (jl.ip is normalized by the condition 

w(0,z)=I m . (1.4) 

The necessary and sufficient conditions for (|1.3p are given in |20(, p. 209] (see also [16 |, 119 ] ). It 
is useful to obtain explicit formulas for H and R. 

The problem (jl.2p is of interest in the random matrix theory: the Markov parameters ap- 
pearing in the series representation w(l, z) = I m + z~ 1 M\ (I) + z~ 2 Mi (£) + ••• are essential for the 
random matrices problems [HE]- In particular, in the bulk scaling limit of the Gaussian unitary 
ensemble of Hermitian matrices the probability that an interval of length I contains no eigenva- 
lues is given by the function P(l), which satisfies the equality A log P(l) = i{rri22(l) — m n(0); 
where mn and 17122 are the corresponding entries of the 2x2 matrix M\. Notice that M\(l) = 
Jq JH(x)dx. 

When J = I m , system (11.2p is essential in the prediction theory |22j . 

We construct a Backlund-Darboux transformation for system (jl.lj) in the next Section 2. 
Section 3 is dedicated to explicit solutions, and Section 4 is dedicated to an inverse problem. 



2 Backlund-Darboux transformation 

To construct Backlund-Darboux transformation we use the methods developed in [14^ [T5] for 
non-isospectral problems and canonical system, respectively. For this purpose fix integer n > 
and n x n parameter matrices A(0), 5(0) = S(0)*. Fix also n x m parameter matrix 11(0) so 



that the matrix identity 

A(o)s(o) - s(o)A(o)* = in(o)jn(o)* (2.1) 

holds. Introduce now matrix functions A(x), S(x) and n(x) by their values at x = and 
equations 

A x = A 2 , U x = -iAUJH, (2.2) 

S x = UJHJ*U* - {AS + SA*). (2.3) 

Then it can be checked by direct differentiation that the matrix identity 

AS - SA* = iUJU* (2.4) 



holds for each x. Notice that the equation A x = A 2 is motivated by the similar equation X x = X 2 
for the spectral parameter A because A can be viewed as a generalized spectral parameter 
(see |14j). In the points of invertibility of S we can introduce a transfer matrix function in the 
Lev Sakhnovich form [HI US US] 

w A (x, z)=I m - ijn(x)*5(x)" 1 (^ - XI n y l Ii{x). (2.5) 

This transfer matrix function has an important J-property |17j : 

wa(x,z)*Jwa{x,z) = J, i.e., wa(x, z)" 1 = Jwa(x,z)* J. (2-6) 
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Put 

v(x, z) = wq{x)~ 1 w a {x, z), (2.7) 
where matrix function wq is defined by the relations 

^wo(x) = Go{x)w {x), w {0)*Jwo{0) = J (2.8) 

G = - J(in*5 _1 n - h jn*5 _1 n + n*5 _1 nj*fl") (2.9) 

up to J-unitary initial value wo(0). (We omit sometimes argument x in the formulas for brevity.) 

Theorem 1. Suppose w is the fundamental solution of system and relations (|2.1|) - (|2.3|) 

are valid. Then in the points of invertibility of S(x) the matrix function 

w(x, z) = v(x, z)w(x, z) (2-10) 

is well defined and satisfies the transformed system 

4-w = iXJHw, (2.11) 
dx 

where 

H{x) = w {x)*H(x)w {x). (2.12) 

Moreover, if det S(x) ^ (0 < x < I) then the fundamental solution of system (|2.1ip is given 
by the formula 

w(x, z) = v(x, z)w(x, z)v(0, z)~ l , < x < I. (2-13) 
Proof. The proof is based on the equation for the transfer matrix function 

-^-wa(x, z) = G(x, z)wa{x, z) — i\wA(x, z)JH{x), (2-14) 
dx 

where 

G(x, z) = i\JH - j(in*5 -1 n - hju*s~ 1 u + n*5 _1 nj*fr). (2.15) 

To prove ([2TT31) consider first ^JITS* 1 . By §F2§ we have 

^JITS" 1 = iJHJU*A*S" 1 - JITS' 1 SxS' 1 . (2.16) 
dx 

Use now ((231) and <t2J6]) to get 

— JITS" 1 = {iJHJ + J)U*A*S- 1 + JITS" 1 ^ - JIT S^UJH JITS- 1 . (2.17) 

dx 

Rewrite identity ((231) as A*S~ 1 = S" 1 A-iS^UJU* S" 1 and substitute this equality into (l2~T71) 
to obtain 

^JITS" 1 = (iJHJ + 2J)ITS- 1 A 
dx 

+ j((hj - i/ m )n*5~ 1 nj - ws^ujhj^s- 1 . (2.18) 
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We shall apply (|2.18p as well as the second relation in (|2.2p to differentiate wa(x, z): 
-^w A = -i(jiJHJ + 2J)n*5" 1 ((^ - XI n ) + XI n ) (A - XQ^U 

+ j({hj - i/ m )n*5~ 1 nj - ics^ujh j)n*5 _1 (^ - Mn) -1 n) 

+ iJIL*S-\A - XI n )-\A x - \ x I n ){A - XlnY 1 !! 

- im*S~ l {A - XI n )~ 1 { - i({A - XI n ) + XI n )UJH). (2.19) 
Use substitutions (A — XIn)" 1 (A — XI n ) = I n and 

(A - XI n )- l (A x - XJ n ){A - XIn)" 1 =I n + 2X(A - XIn)" 1 , 

and collect terms to rewrite (|2. 19|) in the form (|2.14p . 

According to formulas (f277j) (f^9]) . (pJ4"|) . and (j27T5|) we have 

—v(x, z) = wo(x)" 1 (G(x, z) — Go(x))wa(x, z) — iXv(x, z)JH(x) 
ax 

= iXwQ(x)~ 1 J H ~(x)wo(x)v(x , z) — iXv(x, z)JH(x). (2.20) 

Taking into account (|2,8p we get 

w (x)*Jw (x) = w (0)*Jw (0) = J. (2.21) 
Thus we rewrite (|2.20p as 

—v(x, z) = iXJH(x)v(x, z) — iXv(x, z)JH(x), (2.22) 
ax 

where H is given by (pT2|) . From (fl~T|) and (p^2|) it follows that (f2TT]) is true for w of the 
form (I2.10p . In view of (|1.4[) one can see that normalization (|2. 13[) yields w(0, z) = I m . ■ 

Our next proposition provides conditions for invertibility of S. 

Proposition 1. Suppose matrix functions H(x) > and A{x) are summable on the interval 
[0, I], and S(0) > 0. Then S(x) > for < x < I, and so S(x) is invertible. 

Proof. Put 

Q{x) = V(x)S{x)V(x)*, where V x = VA, V(0) = I n . (2.23) 
Then in view of (|2.3p and (|2.23p we have 

Q x = V{S X + AS + SA*)V* = VUJHJ*U*V* > 0, Q(0) = S{0). 
It follows that 

Q(x)>0, S(x) = V{x)- 1 Q(x)(V{x)*y 1 >0. m 

In view of the first equality in (|2.2p invertible matrix function A is of the form A = (B — 
xln)" 1 . Further we shall suppose that A is defined and both A and S are invertible on some 
interval [0, 1]. 
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Remark 1. Suppose A(x) and S(x) are invertible on the interval [0,/]. Using (|2.17p we can 
differentiate 

w A (x, oo) := I m ~ iJ^S^A^U. 

In this way similarly to (|2,14p we can show that the matrix function wq, which satisfies (|2.8p . 
(|2.9p and initial condition ^o(O) = U, admits representation 

w (x) =w A (x,oo)U, U = w A {0,oo)- 1 U, 0<x<l. (2.24) 

Notice that by ([2TT2]) and (pE2"T]) the equality HJH = yields HJH = 0, i.e., if JH is nilpotent, 
then JH is nilpotent too. 

Let now w(l,z) satisfy Riemann-Hilbert equation (|1.2p where W±(s) = lim w(l,s ± irj). 

r/~ »+0 

Suppose all conditions of Theorem Q] are fulfilled. Then, putting 
W±(s) = lim w(l,s±irj), 

we get 

W±(s) = v(l,s)W ± (s)v(0,s)' 1 . (2.25) 
In view of (II. 2|) and (I2.25P we obtain 

W+(s) = W-(s)v(0,s)R(s) 2 v{0,s)- 1 = W-{s)R(s) 2 , 

where 

R(s) = v(0, s)R(s)v(0, s) -1 . (2.26) 

The subcase of nilpotent matrix function R(s) — I m is important [51}]. According to (|2.26p we 
have 

fl(a) - I m = v(0,s)(R(s) - l m )v(0, s)- 1 . (2.27) 
Hence, we get a corollary. 

Corollary 1. If R(s) — I m is nilpotent, then R(s) — I m is nilpotent too. 

3 Explicit solutions 

If we know A, S, IT, then using the results of the previous section we can construct explicit 
expressions for H and R. Consider the simplest case 



m = 2, II S'S. /3= [1 i], 0<x<l, J 



1 

1 



(3-1) 



Then in formula (jl.2p we have 

R(s) = I 2 + vr J(3*(3, 0<x<l. (3.2) 
Indeed, in view of (jl.ip and (|3.ip we get 

f3w x {x, z) = 0, (3Jw x (x, z) = 2i(z — x)~ 1 f3w(x, z), 



i.e., 
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(3w(x, z) = (3w(0, z) = /?, pjw(x, z) = — 2i(ln (z — x))@ + const, 



(3.3) 



where const means some constant (vector). In the first relation in (|3.3p we use normalization 
condition (|1.4[) . Taking into account (jl.4p again, from the second relation in (|3.3|) we derive 



(3Jw(x, z) = 2i ( In 



2 — X 



/3 + /3J. 

Equalities (pOj) and {33]) imply that 

@W+(s) =@W-{s) = /3, PJW + (s) = 2i ( 

+ in ) (3 + /3 J. 



s-Z 



X7T J /S + /9 J, 



2i ( Mb 


s 




s-Z 





From (13. 5j) it follows that 

m + (s) -TVM_(s) = 

Notice that 

TJT* = 2J. 
So according to (13.6H we have 

^ + ( S )-W_( S ) = ^JTV 








" /? * 


4vr/3 




(3J 





47T/3 



2nJ(3*p. 



Moreover, formula (|3.5j) implies that 



W_(s) = -JT*J 
v y 2 



/9 



2i f fin 


s 




s - I 





+ fvr ) (3 + (3 J 



Hence, we obtain 

J (3* = W-(s)jp*. 

Substitute (|3.9|) into (|3.8|) to see that R 2 = I2 + 2i:J[3* (3, i.e., we can assume 
Also we can set 

A = (B - xl n y\ B = diag{6i, 62, • ■ ■ , K}- 
From n x = -iAILJH and ([3TT]> we get 

II J/3* = g = {g k }t =l = const. 
We also have 

-^IL(x)(3* = -2i(B - xI n )- 1 U(x)J(3*. 
ax 

It follows that 

U(x)p* = 2{{ig k ln(6 fc - x)}£ =1 + Zi), h= const. 



(3.4) 



(3.5) 



(3.6) 



(3.7) 



(3i 



(3. 



(3.10) 



(3.11) 
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Formulas (|3.10p and (|3.1ip give us II. We shall assume that b k [0, oo), and so II is well-defined 
on [0, oo). Taking into account (|3.7|) we also get 

U(x)JU(x)* = -IL(x)T*JTU(x)* 

= {{ig k \n{b k -x)}l =l +h)g* +g{{ig k \n{b k -x)} n k=1 +hy. (3.12) 

The matrix function S is easily derived from the identity AS — SA* = zIIJII*. 
Finally, in view of (f!T2lj) and (pTTOj) we get 

0w o (x) = ((3 - ig*S(xy 1 (B - xI n )U(x))U, 
which, taking into account (|2.12p and (13. ip . implies 

H(x) =U*(p- ig*S{x)- 1 (B - xI n )U(x))* ((3 - ig*S(x)- 1 (B - xI n )U(x))U. (3.13) 
From (I2.10p it follows that v(x, s) is J-unitary and we rewrite (|2.27p as 

R(s) = h + v(0,s)(R(s)-h)Jv(0,s)*J. 

Now by (f2"3]l . (|2TTj) . (133]) . and (f^TQ]) we get 

R(s) =I 2 + 7tJw {0)* (f3* - isn(0)*S(0)~ 1 J B(s/ n - B^g) 

x (P* - isUiOyS^BisIn - B)- l g)*w o (0). (3.14) 

Thus matrix functions H(x) and R(s) are given by formulas f)3. 13|) and (|3.14p . respectively. 

Example 1. Consider the simplest case n = 1. Put b\ = b and assume b M. Rewrite (13.121) 
as 



n(x) JU(x)* = i\g\ 2 ( ln(6 - x) - ln(6 - x)) + Kg + gh. 
Here g is the complex number conjugated to g. Hence, in view of (|2.4p we get 

7-n (b — x)(b — x) 



S(x) = (|<7| 2 ( ln(6 — x) — ln(6 — x)) — ihg — igh) 



b-b 



Put now h = to derive 
2i 



S(x) 



b-b 
Rewrite ([37T^|) as 

R( s ) =I 2 + nJU*r(s)*r(s)U, 



= |c?| 2 ( arg(fo — x)) (b — x)(b — x) / 0. 



(3.15) 



(3.16) 



where U = w (0), r(s) = (3 + i bg(S(0)) 1 s(s - 6)~ 1 n(0). By §31$, (IBTTUl) . and (I3TTT) we have 

(3.17) 



n(x) = -U(x)T* JTJ = -g + 2 ln(6 - x)) 1-2 In (6 



2 x ' 2 
Formulas (|3.15p and (|3.17|) imply 

i \ n b—b s 
r(s) =13 + 



46 arg b s — b - 



1 + 2 In b) 1 - 2 In b 



(3.18) 
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Equalities (|3.16p and (|3.18p define R(s) determined by the parameters 6 and g and J-unitary 
matrix U. According to (|2,6|) . (|2.24|) . and (|3.13|) the corresponding matrix function H{x) is of 
the form 

H{x) = U*h(x)*h(x)U, 

where h(x) = [5 — ig{b — x)S(x)~ 1 U(x), and 



U = Jw A (0,oo)*JU = I 2 + 



ib 
W) 



jn(o)*n(o) u. 



Finally, using (|3. 15j) and (|3.17p rewrite h in the explicit form: 

6-6 



h(x) = (3 



l + 21n(6-x)) l-21n(6 



4(6 — x) arg (6 — x) 

4 Inverse problem: explicit solutions 

In view of (|2.6j) it is immediate that formula (|3.14p can be written in the form f)3. 16|) : R(s) 
I2 + irJU*r(s)*r(s)U, where vector function 



r(s) = fiJw A {0, s)*J = fi + isg*(sl n - B*)~ 1 B*S(0)~ 1 T1(0) 
is rational and satisfies the following properties 

r(s) = [n(s) r 2 (s)]eC 2 , r(s)Jr(s)* = 0, r(0) = 0. 
Here J is defined in (|3.ip . Introduce matrices K and j: 



K 



1 

71 



" 1 


-1 " 




' 1 


" 


1 


1 


3 ■= 





-1 



K* = K~ 



KjK* = J. 



Consider function 

u(s) = {V^^ + ^F^)(^(F^-V^F I )y 1 . 

From (|4.2p and (|4.3p we get rKjK*r* = 0, and so 

\u\ = 1, u(oo) = (1 - + 
Rational function u satisfying (|4.5p admits [T7] a so called minimal realization 

u(s) = c 2 (l + ie*SQ 1 {sI n -a)- 1 e), c= (I-O/a/2, 

where 



(4.1) 



(4.2) 



(4.3) 



(4.4) 



(4.5) 



(4.6) 



(4.7) 



Using (|4.6p one recovers H (explicitly, though not necessarily uniquely) from the given function u. 

Theorem 2. Let J-unitary matrix U and rational function u satisfying (|4.5|) be given. Consider 
realization (14.61) and choose vectors 0\ and 62 so that 



:0 1 + c9 2 = O, det (a -icOe^So 1 ) / 0. 



(4- 
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Supposing relations (|4,8p are valid, put 

A(0) :=a-ic69?,S(0)-\ 5(0) := S Q , 11(0) := AK*, A := [9\ 9 2 }. (4.9) 

Introduce now matrix-functions S and U by (|4.9p and equations 

S x = gg* -(AS + SA*), Il(x)Jf3* = g, U x (x)p* = -2i(B - xl^g, (4.10) 

B = A(0)~ 1 , A(a;) = (B-xI n y\ g = 6. (4.11) 

T/ien on the intervals [0, Z], where &et{B — xl n ) ^ and det S(x) ^ 0, matrix function R is well- 
defined by (|3.16|) , (|4.ip . (|4.9|) . and (|4.1ip . Moreover, equality (|4.4|) is true and H corresponding 
to R is given by (|3.13p . (|4.9pH|4,lip , and i/ie second relation in (j2.24|) . 

Proof. First notice that equations (|4.7p and (|4.9p and the first relation in (|4,8p imply identities 
A(0)5(0) - 5(0)4(0)* = iAjA* = TI(0)jn(0)*. 

The correspondence between R and -ff follows now from the results of Section [3j It remains to 
prove ([P]) . For this purpose notice that by (f23J) . ([4TT]) . P~3j) . and (gSD we have 



jiTV^- 1 )* = W(s)jiT/r = W(s) 



(4.12) 



where 2x2 matrix function W is of the form 

W(s) = {W kj (s)}l j=1 =I 2 - ijA*S(oy 1 (A(0) - S I re ) _1 A. (4.13) 
In view of (|4.12p we get 



riis' 1 ) + r 2 (s~ 1 ) cWn(s) + cW 12 (s) 



(4.14) 



ri(s _1 ) -r 2 (s- 1 ) cW 2 i{s) + cW 22 {s) 
From the first relation in (|4.8f) and (|4.13|) follow the representations 

cWu(s) + cW l2 {s) = c - ie* l S{Q)~ 1 (A{Q) - sl n ) _1 0, (4.15) 
cW 21 (s) + cW 22 {s) = c + ie^S(0)- 1 (A(0)-sI n y 1 9. (4.16) 

Using system theory results on the realization of the inverse matrix function, from the first 
relation in (|4.9p and (|4.16p we derive 

{cW 21 {s) +cW 22 {s))~ 1 = c(l -ice* 2 S{Q)- l (a - sQ^O). (4.17) 
Finally, taking into account that i(A(0) - a) = c6 , 6'|5(0) _1 from (|4.15p and (|4.17p we get 

cW n (s) + cWi 2 (s) 2( -0*0-1/ t (A ra x 

cw 21 (s)+cw 22 {s) =c ( i+ies » {sin ~ a) *)■ (4 - i8) 

Equalities (gSD, dH) , and (14TT5D imply ■ 

Compare the result with the recovery of the so called pseudo-exponential potentials (see [5j [7] 
and references therein). 

Remark 2. As |n| = 1, so matrix a in the realization (|4.6p is invertible. Therefore we can 
choose #2 satisfying the first relation in (|4.8|) and sufficiently small for a — icOd^S^ 1 to be 
invertible too. 
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5 Summary 

The first new result in this paper is the construction of the Backlund-Darboux transformation 
for the non-isospectral canonical system (jl.ip . which is important both in prediction theory 
and random matrices theory. The GBDT-version of the Backlund-Darboux transformation, 
constructed in Theorem Q3 is more general than iterated BDT and admits parameter matrix A 
with an arbitrary Jordan structure. (For the applications of GBDT to non-isospectral integrable 
systems see |14].) 

In Section 3, we apply GBDT to the initial system with H = const to obtain a family of 
explicit solutions of system (|1.1[) and of the corresponding Riemann-Hilbert problem (|1.2j) . In 
particular, we construct the transformed Hamiltonians H and the transformed jump functions R 2 
(see formula (|3.13p for H and formula f|3. 14|) for R). The subcase from Example 1 is treated 
in greater detail. The interesting case of non-diagonal matrix A and applications to prediction 
theory will follow elsewhere. 

Finally, in Section 4, using the methods of system theory, we recover H and R from a partial 
information on R similar to the way, in which the Dirac system is recovered explicitly from its 
Weyl function in [7]. 
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